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We compute the model-independent contributions to the primordial bispectrum and trispectrum in 
de Sitter space due to high-energy physics. We do this by coupling a light inflaton to an auxiliary 
heavy field, and then evaluating correlation functions in the Schwinger-Keldysh "in-in" formalism. 
The high-energy physics produces corrections parametrized by H/M where H is the scale of inflation 
and M is the mass of the heavy field. The bispectrum peaks near the elongated shape, but otherwise 
contains no features. The trispectrum receives no corrections at order H/M. 
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INTRODUCTION 

The leading paradigm for understanding the first mo- 
ments of the Universe is known as Inflationary Theory, 
wherein space rapidly expanded in a very brief period of 
time An important property of inflation is that 

its low-energy observables are sensitive to high-energy 
physics [11 [6] . This non-decoupling of energy scales, while 
sometimes confounding, can also be used as a probe of 
ultra high-energy physics such as the effects of quantum 
gravity. In order to verify, constrain or eliminate classes 
of such theories, one must understand the precise rela- 
tionship between low- and high-energy physics in infla- 
tionary theory. 

The primary low-energy inflationary observables are 
the correlation functions of primordial perturbations of 
the quantum field responsible for the inflation. Since such 
perturbations seed the observed temperature fluctuations 
in the cosmic microwave background (CMB), one could 
use the CMB to gain insight into high-energy physics. 
This is particularly exciting given the recent explosion in 
precision data by WMAP \J , Planck Euclid j9j and 
possibly CMBPol/Inflation Probe [TU] . 

The inflationary fluctuation spectrum has been ob- 
served to be nearly Gaussian |7j, and the two-point cor- 
relation probes the free theory: 

(ip^^p) ^ PUk){27rfS\k + p), 

Higher-point correlations then probe the interactions by 
measuring the non-Gaussianity of the statistics. The 
three-point correlation is commonly referred to as the 
bispectrum and is simply the 3-field analog of the power 
spectrum: 

(<y5ki¥'k,<^k3) - B^(ki,k2, k3)(27r)3^3 (j^^ + ka + kg) . 

(1) 

This diagnostic has quickly become an essential tool in 
constraining inflation models [11]. While as yet there 
have only been hints of non-Gaussianity [12l [13] , a defini- 
tive measurement will occur soon from Planck [S]. At 



the most elementary level, non-Gaussianity can be cat- 
egorized according to the shape of the triangle formed 
by the fields' momentum vectors at which the amplitude 
is maximized. The 'squeezed' shape, ki ~ fc2 and 
^3 ~ 0, dominates for local interactions such as <p{x)^ 
|14j . and occur well outside the horizon. The 'equilat- 
eral' shape fci ~ A:2 ~ ^3 occur primarily at the time 
of horizon-crossing and dominates for higher-derivative 
interactions such as found in the Dirac-Born-Infeld the- 
ory [TSKin]. The 'elongated' shape, where ki + k2 ^ fca, 
dominates for an initial vacuum state modified away from 
Bunch-Davies, and often enhances the other types as well 
[T7H^ . A nice heuristic way of understanding this rela- 
tionship is that the more squeezed the triangle, the later 
it originated. In anticipation of the upcoming Planck 
data, a more sophisticated technique has recently been 
pioneered by Fergusson, Liguori, and Shellard [23]. This 
analysis is model-independent in the sense that it does 
not assume a particular template shape (squeezed, equi- 
lateral, or elongated) and so requires an equally sophis- 
ticated model-independent theoretical basis for under- 
standing. 

In order to meet this need, the fundamental question 
how to construct the low-energy effective action in infla- 
tionary backgrounds had to be answered. We did so in 
[21H^ . This new method allows one to reliably compute 
universal corrections to primordial correlation functions. 
The approach is model-independent in the sense that we 
begin with the general form of an ultraviolet-complete 
action and obtain an answer in terms of a small num- 
ber of free parameters that in principle encodes all the 
possible unknown physics. The fundamental organizing 
principle is not yet known. Nevertheless we can simply 
compute the effects in any specific model and general- 
ize from there. We have previously used this technique 
to compute the primordial power spectrum, here we will 
apply it to compute the 3- and 4-point correlations. An 
important difference with flat space effective field theory 
is that a heavy field of mass M will correct the correlation 
function already at order H/M, where H is the scale of 
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inflation. Moreover, as wifl shall show, the leading con- 
tribution to the bispectrum is of the elongated shape just 
as previous studies surmised [T7Hl9l [2T| [22] . 

This article is structured as follows. In §2 we briefly 
summarize the technical framework by which high-energy 
physics can consistently be integrated out of a theory in 
an inflating background. In §3 we calculate the bispec- 
trum corrections, and in §4 those for the trispectrum. In 
§5 we conclude and offer some remarks about detectabil- 
ity and future directions. 



SETUP 
The Generic Action 

We will consider the most general field theory con- 
taining a light scalar field fluctuation if> renormalizably 
coupled to a heavy field %: 
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This is similar to the theory considered in [24j[25], but we 
have added a ip^x interaction to generate non-Gaussian 
interactions, and will now include self-interactions in our 
result [33]. We will take the metric background to be 
exact de Sitter space, 

ds^ = -dt'^ +a{tfdx^, 

where d/a = H — constant. 

The In-In Formalism and Schwinger-Keldysh Basis 

In the Schwinger-Keldysh "in-in" formalism, at some 
early time tin we begin with a pure state |in(tin)), then 
evolve the system for the bra- and ket-state separately 
until some late time t, when we evaluate the expectation 
value 12.7': 



m)) ^ (in(t)|O(t)|in(0) 



(3) 



= (in(ii„)|e'-^4'^*'^(*')0(t)e-'-^4''*"«(*")|in(ti„)). 

Traditionally, the in-state |in) is taken to be the Bunch- 
Davies vacuum state [28], but this is not necessarily so. 
Expanding cosmological backgrounds allow for a more 
general class of vacua, which can be heuristically consid- 
ered to be excited states of inflaton fluctuations. In the 
present context, we will find that integrating out high- 
energy physics generically results in corrections to the 
density matrix, which can be interpreted this way. 

If we denote the fields representing the "evolving" ket 
to be {if+,x+} those for the "devolving" bra to be 



the in-in expectation value ([s]) can be com- 
puted from the action 

S^S[^+,X+\-S[v-.X-]- (4) 

together with the constraint that = It is 

then helpful to transform into the Keldysh basis, 

(p EE + (^_)/2, $=(^+-(^_, 

X = (x++X-)/2, X = x+-X-- 

In this basis the action Q equals 

S[(p, X, X] = - y" d^^dt a{t f [d>fd^ + dxdX + M^^X 
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Fluctuation solutions are easiest written in the conformal 
time T, 

dt = a{T)dT. 

The free field solutions are then 
H 



Vk(r) 



/2fc3 



(1 — ikT)e 



exp 



air) 



1/4 



where for the latter we have used the WKB approxima- 
tion which is always valid for H/M <^ 1. 

We will also need the Green's and Wightman func- 
tions in the Keldysh basis; see [2^ for an explanation 
of their interpretation. Fourier transforming into comov- 
ing momentum, the retarded Green's function can be 
written in terms of the fluctuation solutions, 

G«(Ti,r2) ^ *(^k(ri)$_k(r2)) (5) 

= -20(ri-r2)Im[f/k(Ti)f/^(T2)]. 

The advanced Green's function is then simply the 
time-reversal of this: 

Gi^(Tl,T2) = G^(T2,ri). 

A similar procedure applies for x using its corresponding 
retarded Green's function Q^, 

g^(ri,T2) = z(Xk(Tl)X_k(T2)) (6) 
= ^2ein - T2)lm [Vk(Tl)Vk*(T2)] . 

The Wightman functions are 

^k(Tl,T2) = ('Pk(Tl)'P-k(T2)) 

= Re[C/k(Tl)t/^(T2)], 

= ($k(ri)$-k(T2)), 

•^k(Tl,T-2) = (Xk(Tl)X-k(T2)) 

- Re[Vu{n)VC{T2)], 
= (Xk(ri)X_k(r2)). 
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We will now use these to compute higher-point correla- 
tion functions in the interacting theory. 



BISPECTRUM CORRECTIONS 
Definitions 

A more precise definition of the bispectrum given by 
eq. is 



B^(ki, k2, k3)(2^)^<53 (ki -f k2 + kg) = 

(in(0)|^k,(0)<^k.(0)^U3(0)|in(0)), 
where the time evolution of the in-state is given by 

Mr)) 



(7) 



e J-.n '|m(Ti„)), 



recalling that t — > 0^ corresponds to future infinity. 
From henceforth we will assume a Bunch-Davies in-state. 
There are two types of contributions to 3^^: those arising 
from self-interactions of the </? field, and those mediated 
by the heavy field x coupled to t^. Since they have qual- 
itatively different behavior, we will consider each sepa- 
rately. 



Self Interactions 

Let us first consider the simplest possible interaction: 
that arising from the ip^ interaction. In the Schwinger- 
Keldysh basis this produces two diagrams, shown in Fig- 
ure [ij These are given by (we henceforth omit the 
momentum-conserving delta-function and just assume 
that ki + k2 + ks 0): 

X (^[-iG^^(0,T)]Fk,(0,T)Fk3(0,T) + permutations 
+ \ [-*G« (0, r)] [^^Gl (0, r)] [-^Gl (0, r)] 

Writing out the Green and Wightman functions in terms 
of the functions U\t.{T)^U^{T), and using the fact that 
C/k(0) = C/*(0) = H/y/W yields 



^sclf 



2(27r)3(2fcifc2fc3)3/2i^ 



" dT 



Uk^ir)U^Ar)Uksir) 



+^ki i''')U]i2 {T)Uii3 {''') + permutations 



c.c. 



The integrals will converge in the far past by slightly 
rotating the path of integration into the complex plane, 
so that Tin — — oo(l -I- ie). For the far future, we choose 
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FIG. 1: Lowest-order contributions to primordial bispectrum 
induced by the self-interaction of the light field. Single solid 
lines indicate contractions of i^, dashed single lines indicate 
those of "I>. 
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FIG. 2: Higher-order bispectrum contribution induced by 
multiple self-interactions of the light field. These can be ne- 
glected in slow-roll inflation due to the strong constraints on 
the couplings. 



a cutoff at /i — > , so that the most singular part of the 
resultant integrals are: 



dT 



— ikr 



C30(l+'i€) 



1-n 



(8) 



Thus the late-time correlation is approximately 



5oM 



This correlation peaks when ki ~ 0, corresponding to 
the squeezed shapes of the momentum triangle. Note 
that this is not quite the same local interaction studied 
by Komatsu and Spergel [14j . which introduce nonlin- 
earities into the late-time gauge-independent parameter 
C, whereas we deal with the field fluctuation if and in- 
tegrate the interactions from the far past. Defining the 
non-linearity parameter ^ B^^^ /P^ ^ k^^, the self 
interaction is seen to have a strong red scaling. This 
may possibly be measurable in the near future [3T] using 
the leverage provided by combining the Planck [8] and 
Euclid [9 datasets. Employing two self interactions will 



produce loop corrections of order B] 



self 



5oAoi? and 



DSclf 

B^ 



9'o 



as shown in Figure 2 but slow-roll assumes 
that the these couplings are smaTl and so these contribu- 
tions are subleading. 

High-Energy Correction A 

We now turn to the corrections generated by high- 
energy physics. To leading order in the heavy field, there 
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FIG. 3: Contributions to primordial non-Gaussianity induced 
by the heavy field. The double lines indicate the heavy field 
components {x,X}, with notation analogous to Figure [l] 



are four such corrections, shown diagramaticaUy in Fig- 
ure |3] Since T2 < ti we need only consider corrections 
containing G^{ti, T2) (as opposed to corrections contain- 
ing G^(ti,T2)). 

The first correction is given by 



(-».9i)(-^-^i) 
(27r)3 



dn a(ri)'' / dT2 a(r2)'' x 



•^q-k3(Ti,T2)[-«Gk3(T2,0)] + pcrmutations^ . 

Writing out the Green and Wightman functions in terms 
of f/'s and V^s, we see there are now several types of 
vertices, depending on the permutation of conjugation. 
The first is 



^i(ki,k2) = dr a 



(r)4C/k,(r)C/k,(r)F!(k,+k.)W/W 

° dr (1 - tkiT) (1 - tk2T) . , , 



2\' + 



Then to leading order in H/M the amplitude is 
A(ki,k2) « 



^/(Tc)e'^*Vlki+k2P«?„+i 



2^/h[k2[2kik2{\ - cos e)f' ^ ^/HM 

fci + fc2 + ■\/2fcifc2(l — COS ( 



ki 



(10) 



The physics of this is clear. This diagram accounts for 
the threshold production/decay of heavy particles at high 
redshift in the early universe. Note that in order to evalu- 
ate /(tc), one should use the step-function appropriately 
"averaged" due to the Gaussian fluctuations: 



1 if T > 0, 
1/2 if T = 0, 
if r < 0. 



(11) 



The second possible vertex is identical to Ai but with 
one U conjugated. This has only imaginary-time saddle- 
point solutions. Since our r-integral is confined to the 
real axis we will never pass over this point in our in- 
tegration, and so this amplitude will be suppressed as 
A2 ~ erf(^) ~ ^e^(*^/^)', allowing us to neglect such 
interactions. Finally we consider Az which has both C/'s 
conjugated and so admits no saddlepoint solutions, and 
thus can also be neglected. Note also that while there 
are actually two solutions (differing in sign), we limit 
ourselves to the solution for which Uc < 0, corresponding 
to the inflating phase (see [25] for details on this). 

We can now perform a similar analysis of four-field 
interactions. The first such vertex is 



Sl(ki,k2,k3 



dTa(r)4C/k,(T)(7k,(r)?7k3(T)r 



(ki+ka+ka) 



{r)f{r) 



1 



dr (1 — ifcir) (1 — ik2T) (1 — ik^r) 



A^klklkl y..„ t3 (|ki+k2 + k3 



X exp 



X exp 



-i{ki + k2)T + i j rfr'y |kl + k; 



M2 



where we have introduced a function /(t) to account for 
any step-functions. As first noted in [23], by rescaling 
u = Ht/M the vertex y^i(ki,k2) admits a stationary 
phase approximation at the energy-conservation moment 



i{ki + k2 + k3)T + i / dT\ |ki + k2 + kaP + 



M2 



Using the same rescaling, the stationary phase is now at 

= -V(^l+^2+fc3)2-|ki+k2+k3|2 



which always exists (that is, it is real). To leading order, 
the amplitude is 



ki+k2 = J|ki +k2|2 



(9) Si(ki,k2,k3 



The solution to this defines the New Physics Hypersur- 
face (NPH), 



_ »77r^/(re)e-'^Vl'-^+'-^+'-^l^"..^+^ 
i^HkJ^ 

X [(fci+fc2+fc3)2-|ki+k2+k3|2]"'^' 



' \/2fci/c2(l — cos^), COS 9 



ki ■ k2 

kik2 



ki+k2 + ks + ^/{ki + k2 + kiY - |ki -f k2 + k3|2 



ki+k2+k3|2 



+ u- +\u-„ 
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A similar vertex, but with the third hght field outgoing 
(we place a bar there to remind us of this), is: 

S2(ki,k2|k3) 



dTa(r)4C/k,(T)[/k,(r)l7^3(r)r(k,+k.+k3)M/(^ 



° dr (1 - iklT) (1 - ik2T) (1 + ik^T) 



and the oscillating component ( 12 ) simplifies to 

jj^Ki2/y/2qk3(l-cose) 

where 8 is the angle between — ka = ki + k2 and q. 

Now we must do the integral over q. As done in ear- 
lier work [231 |5S] , we first use the azimuthal symmetry 
around k.3 to simplify the 3d integral into just two vari- 
ables, 



X exp 



-i{ki + fca - ka)T + i J dr'^lki +k2 +k3|2 



M2 



1 



q^dqd{l — cos 9). 



The stationary phase for this is at 

= -v/(fcl+fc2-fc3)2-|kl+k2+k3|2. 

This is not always real, and one must take care. To 
leading order, the amplitude is 



(27r)3 (27r)2 

We wish to transform this into the (u, i?)-basis given by 

u-^ = -y/2k:iq{l-cose), (13) 

E = Mq\u\ = mJ , ^ 

^' ' V 2fc3(l-cos6') 



S2(kl,k2|k3) 



zV^/(T,)e-^^Vlki+k.+k3P.L+i 



^2 



X [(A^i-f fc2-fc3)2-|ki+k2+k3p] 

ki + k2-k3 + ^{ki + k2- k^Y - |ki + k2 -h kgP 



21-3/4 



These can be easily inverted as 

q — ^ , , 1 — cos 9 — 



E 

M\uy 



M 
2k^iEW\ 



so that the q-integral then transforms as 



B 



|ki+k2-hk3|2-Hwr2 + |ui, 

Based on this, it is easy to see that the amplitude is 
51A1 f d\ 



dh 



1 



duEdE 



^ 8(27r)3 (2fcifc2fc3)3/2 J (27r)3 
yli(q,-k3)S*(k2,q,ki) +C.C. 
-A(q, -k3)B;(k2, q|ki) + ^i(q, -k3)6; (ki, q|k2) -|- c.c. 

-|- permutations . 

But the ki -H- fc2 symmetry will cancel the y82-dependent 
terms, leaving only the 61 -dependent terms. 

Examining the loop integral in the expression above, 
the integrand contains a rapidly oscillating component 
due to the difference in vertex interaction times is 



(27r)3 (27r)2M2A:3 
The integrand transforms as 

7r|w|3Me-*7r'^i2" 



(14) 



^i(q,-k3)St(k2,q,ki) 



8HEy%k^ 



f |ki -f k2| + g + v/(|ki + k2| + g)2 - |ki + k2 \ " 
\ ki+k2+q+ ^/{ki +k2 + 5)2 - |ki +\<i2 + qp / 

(12) 

Given such rapid oscillations of this phase, this could po- 
tentially be evaluated using (another) stationary phase 
approximation. However, one can easily convince them- 
self that there exist no such stationary phase points ex- 
cept in the trivial case of fc3 = |ki + k2| = fci + fc2, 
meaning the phase is identically zero. Thus, for all con- 
figurations except this, the amplitude is suppressed as 
H/M. 

To evaluate the integral, we then expand near this 
point, 

K12 = ki + k2- k3, 



0{k^2)- 
(15) 

Noting that _B is a physical rather than a comoving 
scale, we can now easily place limits on the region of 
integration in terms of the physical cutoff A. Recalling 
that at the moment of interaction the energies of the light 
fields will sum to the energy of the heavy field, and so it is 
this latter quantity that we must place the bound of A on. 
The minimal energy bound comes from the geometrical 
constraint cos 9 > —1. The bounds on the loop energy 
are then: 



M 
4kM 
A + VA2 - A'P 
2Mfc3 



< E <A- fc3|u|M, 

A - VA2 - M2 



< u < 



2Mk^ 



The measure (14) and integrand (15) combine to pro- 



duce an £^-integral which is trivial, 

rA-k3\u\M 

dE = A - kMM - 



4fc3|u| 



M 

4^ 



Now including the factor of u^^ from the measure and 
vertex evaluations ( 15 ) we can well-approximate the re- 
maining integral over u as Gaussian near the peak at 
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High-Energy Correction B 



FIG. 4: Example of an 'elongated' bispectrum momentum- 
vector triangle which can be produced from high-energy 
physics, where fci -I- fc2 ~ fcs. 



Mo ~ -3M/8Afc3, 



A /SAfcg 

3" V 3Af 



exp 



M 
Ak3\u 

2^ 



Afcg 

6 V A/ 



(u - UoY 



O 



Extending the domain of integration to u G (—00,00), 
the integral over u can then be easily performed to yield 
another Gaussian for K12. The exponential will cut off 
anything beyond K12 ~ 2'^k^AH/\/3AI^, by which time 
the phase will have completed about half a cycle. Thus 
we can approximate this as a delta-function, 



The next diagram is very similar. 



B 



B _ (-i.9i)(-i-^i) 







% (0, ri)]Fk, (0, n) [-zG«(ri, r^)] x 



(27r)3 



[-«5q_k3(n,T2)][-iG^^(r2,0)] -f permutations 



We do not need to explicitly evaluate this, however, be- 
cause in the stationary phase approximation this is easily 
seen to be identical to diagram A but with a minus sign 
and the Heaviside function 9{ti — T2). This allows us 
to integrate over only half the fluctuations in ri , T2 and 
hence we get 



High-Energy Corrections C and D 

There are two additional diagrams from high-energy 
interactions, 



du e'^fl-^i^^'TTlTrj ("-"0) _|_ 



c.c. 



25 VAfca 



An important caveat is that corrections to this approxi- 
mation will be of order KH/M"^ ^ necessitating that this be 
a small quantity. This is the same restriction one obtains 
for using the plane wave approximation « ^^^"^'^^ 
for all interactions and so is already implicitly assumed. 
The final answer is then 



B^„ « ^9i>^iH^^^ ,5(A:i+fc2-fc3)+perms. (16) 

27y2(27r)4M4fc2fc2fc3 ^ 1 2 3; i V ; 



The requirement that fci -I- ^2 ~ fcs implies that the 
momentum-vector diagram approximates an elongated 
triangle, as shown in Figure [4] The fact that modifi- 
cations to the initial state produce the elongated type 
of non-Gaussianity was anticipated by |17H19j : here we 
see how this originates from fundamental high-energy 
physics. The factor of (A/M)^ can be absorbed into the 
bare couplings 01,^1. Since gi ^ H, Xi ^ 1, A/M ^ 1, 
the correction (16 1 will then scale as B^ ^ H/M. The 

B^/P'^ is scale- 



nonlinearity parameter scaling as /j^l 



invariant, so n 



NG 



0. 



(27r)3 



-zG^^(0,ri)]Fk,(0,Ti)HG^(n,T2)] 



(27r)3 

^c^-\^3{'^l,T2)F■^^{T2,Q) + pcrmutatious 











(27r)3 



(27r)3 



-zG^^(0,ri)]Fk.,(0,Ti)Fq(ri,T2) 



[-«5q_k3 (n ! T2)]^k3 (t2, 0) + permutations 



This time, the stationary phase approximation implies 
the total cancellation of these terms against each other. 
As first observed in the power spectrum calculation [24) . 
this is a general theme: the leading-order corrections aris- 
ing from the dynamics effective action will cancel, and 
the dominant corrections by power counting arise entirely 
from the density matrix. 



Total Bispectrum 



The complete bispectrum correction is then 



Bt 



The high-energy physics produces a scale-invariant elon- 
gated shape of bispectrum corrections. 
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FIG. 5: Examples of 3-point diagrams which are model- 
dependent, and so not considered here. 



There are a few other diagrams involving high-energy 
interactions shown in Figure [s] but they all involve a lo- 
cal or derivative coupling (strongly constrained by slow- 
roll), as well as a factor of H/M. Thus, these model- 
dependent terms will be subleading to the diagrams com- 
puted here. In non-slow roll models local terms may sig- 
nificantly contribute and in fact correspond precisely to 
the terms discussed in [19j . The vacuum modification 
should also enhance local interactions, as first found in 
[3TJ[52]. Higher-derivative actions (for example, DBI [T5] 
or more generally [30]) favor large equilateral bispectra 
[15], although will likely contain new features due to the 
high-energy physics. 




FIG. 6: Contributions to the trispectrum from self- 
interactions. There are also C((7o) and C(Ao) contributions 
which we can neglect in slow-roll inflation due to the strong 
constraints on the couplings. 




B C 



FIG. 7: Tree-level high-energy contributions to the 4-point 
correlation function. 



TRISPECTRUM CORRECTIONS 

The 4-point correlation, or trispectrum, is defined anal- 
ogously to the bispectrum: 

r^(ki, k2, kg, k4) (2^)3,53 (ki + k2 + k3 + k4) 

= (in(0)|^k,(0)^k.(0)^k3(0)¥'k,(0)|in(0)). (17) 

It was suggested in [25] that highly oscillatory corrections 
to power spectrum may be invisible, yet show up in the 
trispectrum as enhanced variance. This is because the 
dominant contribution to the trispectrum is the power 
spectrum-squared, 

^ (5^(ki-k2)(5^(k3-k4)P^(fci)P;^(fc3)-fpermutations. 

A modification in of order H /AI would then produce 
a change of order 2PipH/M in T^p. Below we will eval- 
uate the corrections to T^p, due to self- and high-energy 
interactions. We will see that this is not true. 



Self Interactions 

The first non-factorizable contribution will be from the 
self-interaction as shown in Figure [6| again using a late- 
time cutoff fi ^ 0~ and the most singular part of the 



exponential integral given by (jsj), 

^ (27r)-^ J-oo(i+it) 

(Pki (0, T)Pk. (0, T)Pk3 (0, r) [-iG£ (0, r)] + permutations 
+Pk, (0, r) [^iGl_ (0, r)] [-iG^ (0, r)] [-zG£ (0, r)] 
-|- permutations^ 

iAo dT 

Uti{'^W\^.2{T)U\s_^{T)U]^^{T) + permutations — c.c. 

~ 24(27rfcifc2fc3A:4)^' 

Such local trispectrum correlation functions clearly peaks 
at ki « 0, and has been studied in |32j . This has a 
non-linearity parameter ^ T^°^^ /P^ ^ k~^ scaling. 
There will also be corrections of order T^p ^ OoH^ and 
Tip ^ ^qH^, which are again presumed negligible from 
slow-roll constraints. 



Tree-Level High-Energy Corrections 

Figure [7] shows the three tree-level high-energy dia- 
grams contributing to the trispectrum corrections. We 
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in cosmology) the effective energy encoded in the inter- 
action time is different for xi and X2- A state which is a 
superposition of two different energy states with respec- 
tive energies Ei , E2 and hfetimes ti , T2 famously shows 
oscillations in its time evolution: 



l(x(<)lx(o)>| 



1 



(l + |a|2)2 



-t/n 



(18) 



This is what underlies e.g. neutrino oscillations. To em- 
phasize that it is the same physics here, recall that 



FIG. 8: Loop high-energy contributions to the 4-point corre- 
lation function. 



first evaluate diagram A, 

2 /-O M 

dTi a{Ti)* J dT2 a{T2f x 
([-^G^^ (0, Ti)]Fk, (0, Ti) J-k,+k. (ri, ra) x 



Viki + k2Y - |ki + k2|2 = M/H\ti2\ ■ 
Defining a physical energy similar to eqn (16), 

E12 - {kl+k2)H\Ti2\, 



(19) 
(20) 



can be rewritten as 



rpA 



A/2 



rj.A ^ (-»gl) 

(27r)3 



M2 



cos 



H'^Ti2TzA 

+ In ( 1 



M {Ei2 + M)t34 
ll (£34 + M)ti2 

M T34 
— In — 

H Ti2 



-iGkg (t2, 0)]Fk4 (t2, 0) -f permutations 



£^12 
Id 



In 1 



-E34 



^2 71^4 



4(167r2fcifc2fc3A:4)3/^ 



^l(ki,k2)^t(k3,k4)-f c.c. 



Approximating ln(l+i?/Af) = E/M+ . . . and transform- 
ing to cosmological time r = —e^^*/H, 



2l"7r2(fci/fc2fc3fc4)2M 
[(fcl + fc2)' - |ki + k2n [ih + k,f |k3 + k4p] 



' M2 

-H^^Tl2T34 



COS 



M(ii2 - hi) + 



£'12 

IT 



-E^34 

IT 



permutations 



(21) 



v/(fcl + fc2)2 - |kl + k2 



-iM/H 



^k3+ki + ^(fc3 + fc4)2- |k3-Hk4|2_ 

-|- permutations. 



+ c.c. 



As in the bispectrum, the oscillating phase corresponds 
to the difference in interaction time between the two ver- 
tices. In this case, however, there is no loop-momentum 
integration to wash the oscillations out. That these 
oscillations are really present is in fact easy to un- 
derstand. Due to the cosmological blueshift the dif- 
ference in interaction times translates into a difference 
in energies. The Schwinger-Keldysh-Feynman diagram 
describes a coherent pair of heavy particles sponta- 
neously created each separately decaying into light par- 
ticles at different moments where the on-shell condi- 
tion is reached. One can therefore effectively think of 
the heavy particle as a superposition of "two" states. 



we recognize the familiar AE interference term with a 
new cosmological At term as well. For the specific case 
ki + k2 ~ k^ + ki the energy is conserved, making the 
vertex interaction time equal, and the oscillations disap- 
pear. This is also scale-invariant, gj^^ ~ T^/P^ ^ 1. By 
combining such high-energy interactions with local ones, 
there would be an enhancement as seen in [33]. 

To see whether these characteristic oscillatory features 
are present in the total signal, we must also compute the 
next two diagrams. They are 



{-m? 



^0 

4 



dTi a{Ti)^ / dT2 a(T2) x 



T. 



\x) 



(Ixi) +0i\x2)), one state xi which decays 



(27r)3 

HG^^(0,ri)] Fk,(0,ri) [-z^^ +k. (n, r2)] i^k3(r2, 0)^^k,(r2, 0) 
-I- permutations, 

a(ri)4 / dT2 a(T2)4 [-*G^^ (0, n)] x 



to xi — ^ v{ki) + 'p{k2), the other to X2 fiks) + fiki). 
Moreover (though energy is not a good quantum number 



(27r)3 

Fk,(0,Ti)Hg«+k.(n,T2)] Hg4(t2,0)] [-~iGlir2,0 
+ permutations. 
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Computing their effects, one easily sees that 

The leading H/M contribution therefore vanishes. 



The bounds are now 

M\u\ [i^r' + fcL ^ (fci + fc2)' 

2(fci +k2 + ki2) 
\u±\ ^ 



< E < A-M(fci + fc2)|u|, 



Loop High-Energy Corrections 

Now turning to the loop corrections in Figure [8] the 
first is 



(fci2 + fci + fca) 

fcl2 ± \ 



fc?2 - (^12 + fci + fc2) 



^ (27r)3 

^ ([-iG^ (0, Ti)]i^k. (0, ri)i^q(Ti, r2) x 
-7^ki+k2+q(n, T2)[-iGk3(T2, 0)]^k4(T2, 0) + pcrms.) 



The energy integral is again trivial, and expanding near 
the extremum at |uo| — M/A{ki +^2 + ^12) the resultant 
integrand can be approximated as Gaussian: 



A- M{ki+k2)\u\ 

M\u\ [\u\-^ + kl2 - {ki + k2)''] 



2(fci + ^2 + /C12) 



A^{ki +k2 + ki2) 



I 



27r3()fci)!c2fc3fc4)3/2 



(27r)3 



i3i(ki,k2,q);Bi(k3,k4,q) + c.c. + perms. 



exp 



2M 

-(- 

2\M 



(fci+fc2 + fci2)2(M-U0)' 



There is a rapidly-varying phase present completely anal- Now performing the integral over u, 
ogous to (12), 



' ki + k2 + q + ^/{ki + k2 + q)^ - |ki-t-k2 +qP' 
k3 + k^ + q+ ^{ks + ki + qY ~ |k3 + k4 + q|2 



iM/H 



-«-frKl2M0- 



A(fci +k2 + ki2) 

-0{ki2). 



2 I ffA(fci+fc2 + fcl2) 



As in the bispectrum loop integration, the integration 
over q will only allow contributions near ki + k2 — k^ + k^. 
Now defining variables appropriate to this case, 

K12 = fci + ^2 - (fcs + ^4), 

= -y/{ki +k2+ qY - |ki + k2 + qP, 
E = Mq\u\. 

The [A/MY can be absorbed into the bare coupling 
Taking e to be the angle between fci2 = ki -I- k2 and q, Xi. Unhke the tree-level trispectrum corrections, there 



M 

The final answer is then 
J, _ TrXlH^A^{ki + k2 + ki2] 



T, 



(27r)4Af2(fcifc2fc3fc4)^fcl 



(5(fci-f fc2— fc3— fc4)+perms. 

(23) 



the (u, E) variables can be inverted to yield 
E 



1 — cos 9 = 



M|m|' 
M\u\ 
2k^2E 



is no H/M dependence in the magnitude of (23). Like 
the other high-energy diagrams calculated, it is scale- 
invariant. 

The next diagram is similar. 



[\u\~^ - {k, + k2Y + kl2] 



rpE 

V 



2 /.O 



Just as in the bispectrum case, the integral transforms 
as 



(27r)3 



dri a(ri)** / dra a{T2Y x 



d3q 



1 



duEdE 

(27r)3 ' (27r)2Af2fci"^y 
The integrand transforms as 

7r|M|4M2e-*7f'=i2'' 



d3q 

(27r)3 



-iG^^ (0, ri)]Fk2 (0, Ti)hzG^(ri, r2)] 



(22) [ - zG^^+k2+q(n, r2 )][-*< (r2, 0)]^^k,(T2, 0) -f- perms. 



Bi(ki,k2,q)Si*(k3,k4,q) 



16HEy/kik2k3k4 



+0(ki2). 



It is easy to see that = ~\T/^ . Using the formula for 
Uc in the i3's, we assume that ki + k2 > k3 + k^ to ensure 

Ui> U2. 

Finally, and cancel, and and cancel. 
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Total Trispectrum 



Taking advantage of the cancellations between several 
diagrams, the total trispectrum is 

= {2T:f5\\<ii-\i2)5\]^^-'ki)P^{ki)P^{k^) 
+permutations + T^'="' + +T^. 

While oscillations arise at order H/M in individual di- 
agrams, their cancellation means that no such feature 
appears in the final result. 



CONCLUSION 

We have computed the generic primordial bi- and 
trispectrum corrections in an inflating background due to 
high energy physics. The dominant physics arises from 
the New Physics Hypersurface where the cosmological 
blueshift can just create an on-shell heavy particle. As 
predicted |17H19| . the bispectrum corrections peak near 
the elongated shape of momentum-triangle; here we de- 
rive this ab initio. The bispectrum otherwise contains no 
extraordinary features. The trispectrum has no leading 
order effect. Even though the interference effect of the ef- 
fective mass-shell-time-difference for a pair of heavy par- 
ticles, similar to neutrino oscillations, this cancels when 
all diagrams are taken into account. 

Like that of the power spectrum, the magnitude of the 
bispectrum high-energy corrections are estimated to be 
of order H/M . For optimistic but not impossible val- 
ues of H/M ^ 0.01 these corrections could possibly be 
measured by Planck and other precision experiments 
in the near future. Also like the power spectrum, the 
leading corrections arising from the dynamical part of 
the effective action cancel and the dominant corrections 
arise entire from the density matrix. 

We emphasize that we have only calculated the model- 
independent corrections. Specific models will no doubt 
produce very rich features and should be studied in more 
detail. 
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